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MIMAC : A MICRO-TPC MATRIX FOR DIRECTIONAL 
DETECTION OF DARK MATTER 
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Abstract. Directional detection of non-baryonic Dark Matter is a promis- 
ing search strategy for discriminating WIMP events from background 
ones. This strategy requires both a measurement of the recoil energy 
down to a few keV and 3D reconstruction of tracks down to a few mm. 
The MIMAC project, based on a micro-TPC matrix, filled with CF4 
and CHF3 is being developed. The first results of a chamber prototype 
of this matrix, on low energy nuclear recoils (^H and ^^F) obtained 
with mono-energetic neutron fields are presented. The discovery po- 
tential of this search strategy is illustrated by a realistic case accessible 
to MIMAC. 
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1 Introduction 

, Directional detection of Dark Matter is based on the fact that the solar system 
moves with respect to the center of our galaxy with a mean velocity of roughly 
220 km/s ( Spergel 1988 ). Taking into account the hypothesis of the existence of 
IL^ a galactic halo of DM formed by WIMPs (Weakly Interacting Particles) with a 
„ negligible rotation velocity, we can expect a privileged direction for the nuclear 
' recoils in our detector, induced by elastic collision with those WIMPs. 

The MIMAC (MIcro-tpc MAtrix of Chambers) detector project tries to detect 
' these elusive events by a double detection: ionization and track, at low gas pres- 
sure with low mass target nuclei (H, ^^F or ^He). In order to have a significant 
cross section we explore the axial, spin dependent, interaction on odd nuclei. The 
very weak correlation between the neutralino-nucleon scalar cross section and the 
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axial one, as it was shown in (|Moulin 2QQ5[) ( Albornoz- Vasquez 2011), makes this 
research complementary to the massive target experiments. 



2 Ml MAC prototype 

The MIMAC prototype consists of one of the chambers of a large matrix detector 
allowing to show the ionization and track measurement performance needed to 
achieve the directional detection strategy. The primary electron- ion pairs produced 
by a nuclear recoil in one chamber of the matrix are detected by driving the 
electrons to the grid of a bulk micromegas (IGiomataris 20 06^ and producing the 
avalanche in a very thin gap (128 or 256/im) adapted to the gas pressure. 
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Fig. 1. Track reconstruction in MIMAC. The anode was read every 25 ns in the first 
version, now every 20 ns, giving the 3D track envelope from the consecutive number of 
images, defining the event. The drift velocity of the primary electrons, needed to recover 
the third dimension length, must be determined by an independent dedicated experiment. 



As schematically shown in figure (TJ the electrons are moving towards the mesh 
in the drift space and are projected on the anode thus allowing to get information 
on X and Y coordinates of the recoil track. To determine the X and Y coordi- 
nates with a 100 jam spatial resolution, a bulk micromegas with a 10 x 10 cm^ 
active area, segmented in pixels with a pitch of 350 /im was used, shown in fig. 
2 (center). In order to reconstruct the third coordinate Z of the recoil track, the 
LPSC has developed a self- triggered electronics able to perform the anode sam- 
pling at a frequency of 50 MHz. This includes a specially developped 64 channel 
chip (IRic her 2010[) in its second version, associated to a DAQ (I Bourrion 2010p . 
The electronic card including 8 chips coupled to the anode is shown in fig. [2] (left). 
In total 512 channels are read every 20 ns. 

The 3D tracks are obtained from consecutive read-outs of the anode, every 
20 ns, defining the event. To get the length and the orientation of the track, an 
independent measurement of the drift velocity is needed. An important feature 
concerns the adaptation of the drift velocity value in the gas mixture used to the 
sampling time of the anode. Pure CF4 is too fast to allow to have enough 20 ns 
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Fig. 2. : (Left): MIMAC electronic card with the 8 chips (64 channels each) reading-out the 
512 channels of the lOcmxlOcm MIMAC prototype. (Center): A small portion of the lOcmxlOcm 
pixelized anode. (Right): The bi-chamber module to be installed at Modane in March 2012. 



slices read-out defining a small recoil track (a few mm). We have determined the 
proper mixture to be used in the MIMAC detector coping with very small recoil 
track lengths at low recoil energies (5- 50 keV). We have experimentally demon- 
strated that with a percentage of 30 to 50% of CHFs, depending on the energy 
range of the recoils, we can slowdown the drift velocity of the electrons in the drift 
space, allowing to have at least 2 or 3 read-outs of the anode to define the direction 
of the track. Fig. 3 shows the Magboltz calculation (Biag i 1999D of the drift veloc- 
ity in CF4 as a function of the percentage of CHF^. Preliminary measurements of 
the drift velocity, using an alpha particle source, have been performed in different 
gases in the MIMAC prototype fitting well with Magboltz simulations. 
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Fig. 3. Left: Electron drift velocity as a function of the electric field over the pressure in the 
chamber calculated with the Magboltz code in pure CF4 , in pure CHF2, and in'^He + 5%C4Hio 
mixture for comparison. Right: The electron drift velocity in a mixed of CF^ and CHFs in 50 
mbar as a function of the percentage of the last one for three different drift electric fields 

To determine the kinetic recoil energy of a recoil event we need to know the 
ionization quenching factor (IQF) of such nuclear recoil in the gas mixture. We 
have developed at the LPSC a dedicated experimental facility to measure such IQF. 
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A description of such facility and a precise assessment of the available ionization 
energy by a measurement of the IQF has been performed by (|Santos et al. 2008]) 
in ^He + 5%C4Hio mixture within the dark matter direct detection energy range 
(between 1 and 50 keV) . The IQF measurents of ^^F in CF4 has been performed 
recently (|Guillaudin 201 1|) giving the recoil energy threshold to be expected in 
directional detection using a ^^F target. We have detected ^^F with 5 keV kinetic 
energy leaving in ionization only 1.2 keV. A publication concerning these important 
measurements for dark matter direct detection, as a funcion of the gas pressure 
and up to 50 keV, is in preparation. 

3 Experimental results 

In order to validate the MIMAC prototype we had to be able to detect simultane- 
ouly the energy and the track of nuclear recoils in the range of 1- 100 keV. Neutrons 
are the ideal tool to produce nuclear recoils by elastic collision as the Wimp event 
signal searched. We have performed dedicated experiments at AMANDE facility 
(|Lebreton 20111) (IRSN - Cadarache (FRANCE)) producing neutron monochro- 
matic fields from (p,n) reactions. The neutron energy is selected from the maxi- 
mum energy produced at the resonance at 0^ by a moving arm around the target. 
In such a way, neutron energies covering 10 to 565 keV are available. Two different 
neutron fields were used at 144 keV and 565 keV maximum energies on C4H1O and 
CF4 respectively. A X-ray calibration has been produced from ^^^Cd and the ^^Fe 
radioactive sources giving the 3.05 keV and 5.9 keV lines respectively. The X-ray 
spectrum is shown in fig. 4 for a high gain run. This spectrum is the projection 
of all the events producing tracks in the active volume. 



Fig. 4. X-ray calibration with ^^^Cd and ^^Fe sources. We can see the 3.05 keV and 
the 5.9 keV lines. The second one broadened by the K/s 6.4 keV not resolved. 
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A very important point for dark matter detection is the ability to separate 
gamma ray background from nuclear recoils. For a given energy, an electron track 
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is an order of magnitude longer than a recoil one. In such a way, the discrimination 
of electrons from nuclei recoils by using both ionization energy and track length 
information, is performed in its first step. The plot of the length of the tracks 
versus their ionization energy, shown in fig. 5 (left) was the first raw analysis during 
the experiment with 144 keV neutrons on C4H10, showing the high quality of 
observables available from the read-out. This discrimination is even improved at 
low energies event by event by a new degree of freedom called the normalized 
integrated straggling (NIS) of the track, defined as the sum of all the angular 
deviations along the track, normalized by the ionization energy of the particle. The 
NIS histogram shown is produced using the regions of the length vs. energy plane 
where no confusion arises between electrons and recoils giving for discrimination a 
clear cut on the NIS value. This NIS discrimination degree of freedom is available 
from linking the barycentres of each 20ns slice read-out, putting in evidence the 
difference in the collision straggling between electrons and recoils. The track of 
recoils are more straight than those of the electrons due to their larger mass, as is 
illustrated in fig 5 (right). 




Fig. 5. : (Left): The length vs. ionization energy of all the raw events produced by the 
neutrons (144 keV) on C4H1O and gamma rays produced on the target. (Right): The NIS, 
obtained as a sum of all the step deviations (every 20 ns) normalized by the energy, in 
black for recoils and in blue for electrons. 



In order to understand the wealth of the information obtained by the read- 
out, we show in fig. 6 the 3D reconstruction of a proton and fiuorine tracks 
showing the shape of each X-Y consecutive read-outs. A straightforward idea of 
the direction of the track is yielded by the barycentre of each patch. The shape 
of each patch gives information on the direction of tthe track inside the 20 ns 
drift, that is extracted by a maximization of a likelihood in the final step of the 
analysis (Billard ~et al. 20T1 1) giving at the same time an estimation of the distance 
from the anode of the track, taking into account the longitudinal and transverse 
difusion of the primary electrons opening the possibility of a fiducialization of the 
active volume. 
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Fig. 6. Left: A 57 keV 3D proton track in CF4 plus 30% CHF3 in 50 mbar. Right: A 
40 keVee fluorine track in the same gas mixture 



4 MIMAC prospective 

The MIMAC first bi-chamber prototype consisting of two chambers of 10x10x25 
cm^, shown in fig. 2 (right), will be installed at the Laboratoire Souterrain de 
Modane (LSM- France) in March 2012. We are working on the 20x20 cm^ pixelized 
micromegas having 1024 channels of read-out. The demonstrator of Im^ should 
follow the validation of the first bi-chamber prototype. The Im^ will be built as the 
reproduction, 100 times, of the bi-chamber (20x20x25 cm^) module, divided in two 
big chambers 0.5 m^ active volume each. This should be follow by a large matrix 
detector 50 m^. We have shown recently that the directional detection opens the 
way for a discovery with a reasonable number of WIMP events (jBillard 2010ap . 
Figure 7 presents, in the proton axial cross-section versus WIMP mass plane, the 
discovery region and the background free exclusion limit for a 50 m^ during 3 
years. The Im^ demonstrator detector may allow for a discovery down to 10 ~^ pb 
and for an exclusion down to 4. 10~^ pb. 
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Fig. 7. Axial cross section on proton (pb) as a function of the WIMP mass. The black 
curve shows the 3a discovery region. The shaded region shows the difference with respect 
to the best exclusion curve today (|Behnke 2011|) in the spin-spin interaction vs. WIMP 
mass plane. The exclusion curves available at such exposure are shown in red. The red 
dashed curve shows the exclusion curve with 10 background events / ( kg year) contami- 
nation. The green region represents the SUSY models calculated by (^Trotta 2011)) . 
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